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Abstract 

We have studied various classical solutions in E? cosmology. Especially we have obtained 
general classical solutions in pure cosmology. Even in the quantum theory, we can solve 
the Wheeler-DeWitt equation in pure cosmology exactly. Comparing these classical and 
quantum solutions in R^ cosmology, we have studied the problem of time in general relativity. 
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1 Introduction 



Recently quantum theoretical treatment of gravity has been paid attention in various 
topics such as the time development of quite early universe |]I], 0], exact treatment of Hawking 
radiation]^ and the problem of time in cosmology 

On the other hand, attempts to generalize Einstein's gravitation theory have a long 
history in various purposes. The generalized Einstein theory with additional term was 
introduced to regularize ultraviolet divergences in the Einstein theoryp]. It was applied 
to cosmology to obtain the bounce universe to avoid the singularity at the creation of the 
universe0. The structure and the properties of the higher derivative gravitation theory were 



further elaborated in subsequent works |]I0 



Other people have tried to generalize Einstein theory to modify Newtonian force laws from 
the framework of general relativity |TT], |T^. Since there is no established direct observation 



of dark matter, there are many attempts to explain rotation curves without dark matter 
by modifying the Newtonian force empirically [|T^] or from the framework of the general 
relativity [0]. Whereas, on the application of the generalized Einstein theory to cosmology, 
many people are interested in cosmology as it may give a promising model to explain the 



inflation in early universe without using the fictitious scalar field|15, IG]. The semiclassical 



solution of higher-order cosmology gives the de Sitter universe jlTl] , but there arises instability 
in the quantum treatment [|18| . 

In the previous paper, we attempted to explain not only the fiat rotation curves of spiral 
galaxies but also the large-scale structure of the universe, starting from a simple model with 
adding R^ term to the Einstein action|l^. For the generalized Einstein action, we have the 



quite stringent observational constraint that the coefficient 7 of the Robertson expansion [P^j 
must be quite near 1 from the classical tests of general relativity by the radar echo delay 
experiment. However, in the previous paper, we constructed the model which gives 7 = 1 
at solar scale and gives Sander's type gravitational potential at galaxy scale pO|. 

In series of our papers, we applied our R^ involved gravitational theory to 1) solar 
scale, 2) galaxy scale, 3) large-scale of the universe, such as "Great Wall", in classical 
level. Then in this paper, we will apply this R^ involved gravitational theory to cosmology. 
First we will find various classical solutions in R^ cosmology. Especially we will give general 
classical solutions in pure R^ cosmology. Next we will consider the quantum theory of pure 
R^ cosmology, and will give exact solutions of Wheeler-DeWitt equation. Comparing these 
classical and quantum solutions in R^ cosmology, we will study the problem of time in the 
general relativity according to the de Broglie-Bohm interpretation [I 
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2 Classical Solutions of Cosmology 



We use Weinberg notation |jT9| and consider the following generalized action: 



/ = J d'^Xy/—g 



16nG 



matter 



where G is the gravitational constant, R is the scalar curvature and a is a constant. 

In order to study classical and quantum solutions in R^ cosmology, we consider the 
homogeneous and the isotropic space-time with the metric 

ds^ = -hiifde + aitf [dr^ + r\de^ + sine^d^^) 

where we put the curvature term to be zero. 

Using this metric, we obtain the expression of the scalar curvature 



R 



( d c? db 

-6 \ 

V ab"^ a'^b'^ ab^ 



(3) 



After partial integration, the original action can be transformed into the form 



d X / dt 



ad^ „ ( 'a^a 2ddba d^b'^a d'^ 



Qa 



+ 



65 + ab^ 



+ a^bC 



matter 



SttG [ b ' "'^\b^ b^ 
For perfect fluid approximation of matter, equations of motion to be solved become 



(4) 



H'^ - 6a [2HH - + 6HH^} 
p + 3H{p + P) = 0, 

where we put 6 = 1 after variation. We use notation H 
density and pressure of perfect fluid respectively. 



8nG 



(5) 
(6) 

d/a, and denote p and P as 
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2.1 Classical solution with perfect fluid in pure R'^ cosmology 



We put P = 7p in this case, then equations to be solved become in the form: 

-6a {2HH -H^ + 6HH^} = ^p, (7) 
p + 3/f(l + 7)p = 0. (8) 

From Eq.(^, we obtain pa^^^^^^ = Poa^'^^^^'''^ = (const.). Substituting this relation into 
Eq.(|3); we obtain the following power type solution of the universe: 

a{t) = kit^^, (9) 

"^'^ = "° [ 4a(5 - 37) ) ' 
where we use the initial condition a(0) = 0. 



2.2 Classical solution without matter in R + aR^ cosmology 

a) Case of a < 

In this case, equation to be solved becomes in the form[^, T^]: 

H'^ -Qa[2HH - + 6HH'^} = 0. (10) 
When we denote a = —a{> 0) and change variables into forms 

H{t) = i^mr, 



t = (11) 

we obtain following equation: 
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(12) 



except the uninteresting trivial case / = (static universe). 

Though we cannot find exact solution for Eq.([12D, we can solve that equation perturba- 
tively. First, treating the third term in Eq.(|T^) as the perturbation, we obtain the following 
power series solution: 

3^ 25^^ 4725^^ 4606878^^ 4558125^ ^ ' ' J ' ^ ' 

where we assumed that singularity occures at ^ = 0. 

Second, treating the second term in Eq.(|l2D as the perturbation, we obtain the following 
oscillating solution by taking the initial condition /(O) = 0: 



/(0 = sine + 



1 ^ 

— (cos ^ — COS 3^) sin ^ 

32 8 



+ 



(14) 



Combining these perturbative solutions, we know the asymptotic behavior f{^) ~ sin(,^ — 
^o)/v^ as ^ ^ oo. 

Here we must notice that the above solution does not approach the solution of equation 
= 0, which comes from Eq.(|l^) by putting a = from the beginning, even in the limit 
a 0. This is because a is the coefficient of the highest derivative term. Then the limit 
(5-^0 give the quite chaotic behavior, which is quite similar to the turbulent limit in the 
Navier-Stokes equation. Actually, H{t) has the following t dependence: 



H{t) = n^). (15) 

36v« v24a 

In the limit a with finite t, H{t) oscillates quite rapidly and the oscillation amplitude 
diverges, and finally becomes chaotic. 



b) Case of a > 
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We define a = a{> 0) in this case and use the same notation as those of Eq.(|T2|), then 
we have the following equation: 



f + 4-/^0. (16) 

We cannot solve exactly in this case either, but we can see that /(^) ~ \/ — ^o) as ,^ — > oo. 



2.3 Classical solution without matter in pure 

In this case we put H{t) = f'^{t), then equation of motion to be solved becomes in the form 
£ f dp 

—— H — — = 0, in addition to the trivial case f = 0. The first intee;ral of this equation eives 

df n o 

— + / = — Ci , where ci is the integration constant. The special solution of this equation is 
dt 

given by / = — Ci = (const.), which gives the de Sitter universe a{t) = a(0)expci^t. Most 
general solution is given in the following implicit form: 

1 , ( {f + c,f \ ^ 1 ^ 2/-C1 



In ci = case, we can express / in the explicit form / = l/y2(t — to), which gives the 
square root type solution of the universe 



a{t) = a^^/t — to. 



3 Another Formulation of Cosmology 

In order to reduce the higher derivative term, we introduce an auxiliary field a and 
rewrite the original action in the form: 
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d Xy 



R + a{-2Ra - a^)] +C 



(19) 



This reproduces the original action after substituting the solution of a field. In the quan- 
tum theory, taking the path integral formulation, we have additional term 1/ [^di'^)^^^ 
when we integrate over a. This additional factor changes the measure of the inner prod- 
uct, which is connected to way how to define the ordering of operators. In this sense, our 
quantum theory of gravity and conventional quantum theory of R^ is different only in 
the measure. As there is no a priori criterion to determine the measure, here we take above 
action as the starting action in quantum gravity. 



Using the metric of Eq.(0), we obtain 



(Tx / dt 



aa n, / 

6— + aa^h I2a 



a a 
+ 



x62 



62 




+ a^hC 



matter 



.(20) 



In pure R^ cosmology, we can analize exactly. Then we will consider only pure R^ case 
hereafter. In pure R^ cosmology, the Lagrangian density becomes 



aa^h ( ( d a} 
-12or| — + 



ah 



a 



(21) 



IGttG V \a}P a262 afes^ 
We define X+ = a^/^a^/^,X_ = a^^'^,B = ho^l'^ , and after partial integral, we obtainp 



a 
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Using this form of Lagrangian, we obtain equations of motion 



(22) 



x+x_ — ~rc~"^+"^-' 

lo 

^± = (23) 

Combining these equations, we can obtain general classical solutions in pure i?2 cosmology 
by the following two steps. 
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step 1) 

First we use temporal time variable T, which is defined by the condition B{T) 
solve equations of motion with this T. The metric at this step gives 



= 1, and 



ds' = -b{TydT' + airy [dr' + r'de' + sin e'd^'\ . (24) 

Step 2) 

The original time variable t and the temporal time variable T is connected with the 
condition dt = b{T)dT, then the classical solution of the original metric is given by a{t) = 
a{T{t)). 

In Step 1, equations we deal with are foUowings: 

dX+dX_ 3 ^ ^ 



dT dT 16 



d^X+ 3 



which come from Eq.(^) after substituting B = 1. Then we get general solutions of the 
following three types: 



solution 1 

X+ = X+(0)expcjT, 

X_ = X_.(0)expcjT, (26) 
solution 2 

X+ = X+(0) coshcj(T - To), 

X_ = X_(0) sinhcj(T - To), (27) 
solution 3 

X+ = X+(0) sinhcj(T-To), 

X_ = X_(0)coshcu(T-To), (28) 



with uj = y3/16. From these quantities, we construct a(T) = (X_)^/^,6(T) = {X^/X^Y^^. 

In Step 2, solutions is expressed with the original time variable t via the temporal time 
T as follows: 
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solution 1 




t-to 




a{T)- 


= ix^ 


solution 2 




t-to 


= {X- 


a{T)- 


= iX- 


solution 3 




t-to 


= {X- 


a{T)- 





(29) 



(30) 



(31) 

— oo reproduce the previous 
solution a{t) = aoVt — to given in Eq.(|TBp. 

General solutions in these forms are more useful than the expression in Eg. {^i7\) . This is 
because f(t) is not given in the explicite form in Eq.(|l7D, and also because we must integrate 

2, da(t)/dt 



In the above solutions, X±{0) — > (finite) x expcuTo with Tq 



further to obtain a{t) from / (t) 



a{t) 



even if we know the explicite form of /. 



4 Quantum solution in cosmology 



When we solve classical solution in pure gravity, we noticed that the auxiliary field 
formalism is quite powerful. Therefore, even in the quantum theory, we expect that it will 
become easy to solve the Wheeler-DeWitt equation when we use auxiliary field formalism. 

Using the auxiliary field formalism, we start from the previous Lagrangian Eq.(|22D 



L 3-D 



f32) 



with a = a/lGirG. We define the canonical momentum in a standard way 
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p_ 



dC 
dC 



16d • 



dX_ 3B 
dC 

Pb = — - ~ (primary constraint). 
dB 



The Hamiltonian is given as follows: 



n = P+X+ + P_X_ - £ 



35 
16d 



-P,P_ - aBX:X. 



The secondary constraint, which comes from Pb ~ 0, gives 



(33) 



(34) 



PB = {n,Pb} 



cl. 



16a 



-P~P^ - aX,X_ ^ 0. 



Taking the Schrodinger representation P+ = —i 
Wheeler-DeWitt equation of the form 



d 



dX, 



P_ = -i- 



d 



dX_ 



EqM 



(35) 
gives the 



+ ——x+x. 



dX+ (9X_ 



vI/(X+,X_) = 0. 



(36) 



In the above equation, the operator in square brackets is symmetry under the transfor- 
mation 



X- - ]x., (37) 

where i is an arbitrary constant. If we find a solution \E'(X+, X_), \E'(£X+, 'j-XJ) also becomes 
a solution by using the transformation of Eq. (|37D in Eq.(^). 
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It is convenient to rewrite with variables Xq, Xi, which are defined by X^ = j{Xq + Xi] 
X_ = £{Xq — Xi). Using these variables, Eg . (|36D becomes in the form 



dXl dXi 



{xl - XI 



xl>(X+,X_) = 0. 



(38) 



In order to find the eigenstate of the state vector, we put the state vector into the sepa- 
rable form \1/(X_|_,X_) = \l/°(Xo)\l/^(Xi). Using this separable form, the above equation be- 
comes separated into two Wheeler-DeWitt equations, which are nothing but the Schrodinger 
equation of the harmonic oscillator with upside-down potential [^. 



Q2 



dxi 



64a2 
3 

64d2 



dXf 



xl] ^\X,) 



E^\Xo) 
E^^(Xi) 



Solutions are given by 



^°(Xo) = ^^(Xo), ^i(Xi) = <l>„(Xi) 



with 



(39) 



(40) 



<I>„(X) = H4e-'^/'(3X)e 



E = Er, 



-2il3^{n + l/2), P 



64d 



2\ 1/4 



(41) 



These $„(X) are nothing but analytic continued forms from those of the harmonic oscillator. 
As equations are linear, superposition are also solutions. General solutions are given in the 
form 



v^(X+,X_) = Y: [c.$.(Xo)$.(XO + ci.$:(Xo)$:(XO] 



(42) 



n=0 



Here we must comment on these solutions. First of all, these state vectors are not normal- 
izable. The origin of this is because of the instability of the gravitational theory. According 
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to the standard Copenhagen interpretation of the quantum theory, the above solutions are 
unphysicaL But the above solutions are realy quite physical when we compare the above 
quantum solutions with classical solutions, as we will discuss on that later. Therefore, here 
in this paper, we take milder standpoint than the standard Copenhagen interpretation es- 
pecially in quantum cosmology, and allow all solutions of the state vector without imposing 
the normalizability. As we do not impose the normalizability on the state vector, we obtain 
twice as many solutions compared to those of the harmonic oscillator. 



4.1 Correspondence between classical and quantum solutions un- 
der de Broglie-Bohm interpretation 

In this section, we study the connection between classical and quantum solutions to un- 
derstand the problem of time in the general relativity. In the general relativity, because 
of the invariance under general coordinate transformation, the consequent Hamiltonian be- 
comes weakly zero, then we come to the problem of time in canonical formalism. This is the 
problem because we may naively think that the weakly zero Hamiltonian will give only the 
solution of unrealistic static universe. 

To overcome this unrealistic conclusion, many people take an approach from the WKB 
approximations [^, |^, 0. The essential idea of the WKB approach is the followings. We 
first separate the Hamiltonian into classical and quantum parts, and we introduce time 
through the classical solution. Then the quantum Hamiltonian plays the role of the time 
development operator in the quantum system. Though classical and quantum Hamiltonians 
are both time developing, total Hamiltonian balanced to becomes zero. If there is matter 
field in the system, this approach may be reasonable, as the classical Hamiltonian may be 
dominated from the metric contribution and the quantum Hamiltonian may be dominated 
from the matter field. We cannot see the time dependence of the metric on the co-moving 
system, the matter may develop with time according to the quantum Hamiltonian. 

We will not take this standard approach in our pure cosmology. Because there is no 
matter, it is unrealistic that classical and quantum contributions from the metric have the 
same magnitude. In the following, we take the de Broglie-Bohm interpretation |]2T[] between 
classical and quantum solutions and try to study the problem of time. 

In order to compare classical and quantum solutions in cosmology, we first consider the 
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simplest state \l/o. The state vector in this case is given by \I/o = e'-'^^^^o'^^^^^'^ = e*^^('^++'^-)/'^. 
In general, we decompose the state vector into radial and phase components in the form 

/?2 1 

* = Ae'^. Then S for n = state is given hj S = -^{fXl + -^Xl). 

According to the de Bloglie-Bohm interpretation, we obtain the first order differential 
equations for X± as 



16a • dS ei3^X^ 

-A_ — P_i 



3 " ^ dX+ 2 ' 
16«^ _p _ dS _P'X. 



Prom the above equations and noticing relations = Sa/v^, = y3/16, we obtain Xj^ = 

i~'^ujX^, X_ = £^ti;X+. Then we have the second order differential equation X± = lJ^X±. 
The general solution is given by 

X_ = ^!^=£(cie-^-C2e-'^^). (44) 

CO 



If we compare these solutions with those of Eq. (|25D , the above solutions correspond to 
the classical solutions, Eqs.(PBp- (ESI), exactly. The more detailed correspondence between the 



parameters in Eq.(44) and Eqs. (|2q) - (|28|) is given as follows with £ = wX_(0)/X+(0): 



solution 1 (ci > 0, C2 = 0) 

ci = v/x+(0)X_(0), 
solution 2{ci > 0, C2 > 0) 

ci = v'X+(0)X_(0)e-^«, C2 = v'X+(0)X_(0)e-^°, 
solution 3{ci > 0,C2 < 0) 

ci = v/x+(0)X_(0)e-"^°, C2 = -7x+(0)X_(0)e"^«. (45) 

Because of instability in cosmology, that is, the potential of the harmonic oscillator 
being upside-down, the solution becomes the sum of exponentially blowing up and damping 
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term as time developes. This behavior of the classical solution is quite physical, we relax 
the standard Copenhagen interpretation and allow the unrenormalizable state vector in the 
quantum cosmology. 

For general state vector with n > 2, expressions obtained from the de Broglie-Bohm 
interpretation is different from expressions of classical solutions. This is because only n = 
and n = 1 state vecotor is composed of only one term. 

If the arguments of the Hermite polynomials become large, which describes that the 
universe becomes large, expressions constructed from the de Broglie-Bohm interpretation 
approach classical solutions. This is because the additional phase of wave functions takes 
constant value —mi/2 in the above case, and constant phase does not contribute in dS/dX±. 

According to de Broglie-Bohm interpretation, time comes through the relation 
X = P = dS/dX, so that even if the Hamiltonian becomes zero system develops accord- 
ing to time. 

In the rest of this section, we will point out the relation between de Broglie-Bohm in- 
terpretation and Ehrenfest theorem^, which describes classical orbits as expectation values 

1 92 

in quantum mechanics. For the one dimensional Hamiltonian H = 2 + the 

2m dx 

Ehrenfest theorem gives 



after partial integration. While, according to the de Broglie-Bohm interpretation, we obtain 



m- 



dX 
It 



p \ dx dx , 



x=X 



dx 



(47) 



x=X 



If we compare Eq.(^) with Eq. (^Tf) , we can understand that de Borglie-Bohm interpretation 
is the local form of Ehrenfest theorem. Here we emphasize that the state vector is assumed 
to be normalizable in the proof of Ehrenfest theorem, while the normalizability of the state 
vector is not assumed in the local form of Eq.(^). In fact, the state vector is not normalizable 
in our examples. 



14 



5 Summary 



In this paper, we have studied various classical solutions in the generalized Einstein cos- 
mology which contains higher derivative terms R^. Especially we obtained general classical 
solutions in pure cosmology. Even in the quantum theory, we can solve the Wheeler- 
DeWitt equation in pure R^ cosmology exactly. We compared classical and quantum so- 
lutions in this R^ cosmology, and studied the problem of time using the de Broglie-Bohm 
interpretation. Further we pointed out the relation between de Broglie-Bohm interpretation 
and the Ehrenfest theorem. 
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